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Generating Entangled Two-Photon States with Coincident Frequencies 
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It is shown that parametric downconversion, with a short-duration pump pulse and a long nonlin- 
ear crystal that is appropriately phase matched, can produce a frequency-entangled biphoton state 
whose individual photons are coincident in frequency. Quantum interference experiments which 
distinguish this state from the familiar time-coincident biphoton state are described. 
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Spontaneous parametric downconversion (SPDC) has 
been the entanglement source of choice for experimental 
demonstrations of quantum teleportation, entanglement- 
based quantum cryptography, and Bcll-incquality viola- 
tions, etc. However, the biphoton state generated via 
SPDC under the customary phase-matching conditions is 
maximally entangled only when a continuous-wave (cw) 
pump is used |l|, §]. In pulsed-pump experiments, the 
fringe visibility in biphoton interference measurements 
decreases as the duration of the pump pulse is reduced 
H . The timing and pump-intensity advantages of pulsed 
experiments has thus spurred work to retain or restore 
maximal entanglement in the pulsed regime [Q. 

In this paper a new method for obtaining maximal 
entanglement from pulsed SPDC is reported. Our ap- 
proach uses a long nonlinear crystal and extended phase- 
matching conditions tailored specifically to pulsed op- 
eration. As such, it does not require any filtering or 
post-selection, and it reaps the high-conversion-efHciency 
advantage that long crystals afford. Furthermore, the 
biphoton states that it produces are comprised of pho- 
tons that are coincident in frequency, in contrast to the 
usual cw phase-matching case whose biphotons exhibit 
coincidence in time. Coincident-in-frequency entangle- 
ment is important because the iV-photon version of such 
a state has been shown to improve the accuracy of time- 
of-flight position sensing or clock synchronization by a 
factor of VN §. 

Consider SPDC with a cw pump and conventional 
phase matching that is operated at frequency degener- 
acy, i.e., phase matched such that the center frequencies 
of its signal and idler equal half the pump frequency. This 
system produces a biphoton of the form, 



ITB) 



2^ 



(w)|w p /2 - uj) s \uj p /2 + uj) t 



(1) 



Here: \u) s ) s and are single-photon signal and idler 
states in which the photons are present at frequencies uj s 
and Wi, respectively; u> p is the pump frequency; and <fi(uj) 
is the spectral function of the state, so that \4>{u) p /2 — ui)\ 2 
is the signal's fluorescence spectrum. The notation |TB) 
indicates that Eq. (Q) is the usual twin-beam state of 



SPDC. The frequency entanglement of this state dictates 
that a signal photon at frequency w p /2— uj is accompanied 
by an idler photon at frequency u p /2+lj. The sum of the 
signal and idler frequencies is therefore fixed at the pump 
frequency. By Fourier duality, this implies that the signal 
and idler photons occur in time coincidence — to within 
a reciprocal fluorescence bandwidth — as has been shown 
in the famous "Mandel dip" experiment Q . 

On the other hand, the SPDC biphoton that will be 
studied in this paper is, 
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In this state, a signal photon at uj p /2 + lu is accompa- 
nied by an idler photon at the same frequency. This 
coincident-in-frequency behavior leads, via Fourier du- 
ality, to symmetrically located occurrences in time. In 
particular, a signal photon appearing at To + t is ac- 
companied by an idler photon at To — t, where To is 
the mean time-of- arrival of the biphoton pulse. Because 
this biphoton possesses a narrow distribution in signal- 
minus-idler difference frequency, we have dubbed it the 
difference-beam (|DB)) state. Note that its mean time- 
of-arrival, To, plays the role of the fluorescence center 
frequency, w p /2, in comparing the DB and TB states. 
Thus, whereas the photons in |TB) may be discriminated 
by frequency measurements, those in |DB) may be dis- 
tinguished via time-of-arrival measurements. 

The rest of the paper is organized as follows. First, we 
derive the output state of SPDC. The phase-matching 
conditions that are needed to create the DB state are 
then obtained and explained. Next, we present quan- 
tum interference experiments that can distinguish be- 
tween the TB and DB states. Finally, we give a feasi- 
bility study for |DB) generation using periodically-poled 
potassium titanyl phosphate (PPKTP). 

A textbook treatment of the SPDC process (see for ex- 
ample Q]) allows us to deduce the state at the output of a 
compensated SPDC crystal. We will give a brief deriva- 
tion here assuming colinear plane-wave operation. In the 
interaction picture under the rotating-wave approxima- 
tion, the Hamiltonian that gives rise to the creation of 
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the two downconverted photons starting from a single 
pump photon is given by, 



Hi(t) = 

j-L/2 

S / dz x (2) 4 +) ( z > *)^ _) ( z > ( z , t) + h.c 

J-L/2 



(3) 



where 

X (2) is the nonlinear coefficient and L is the length 
of the downconversion crystal, S is the pump-beam area, 
and and E^> = (E^Y are positive-frequency and 
negative-frequency electric field operators with the sub- 
scripts {p, s, i} denoting pump, signal, and idler, respec- 
tively. These electric field operators obey, 
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for j = p, s, i, where aj(u)) is the annihilation operator for 
frequency-w photons, Tlj(u>) is the refractive index for the 
jth beam (pump, signal, or idler), and kj(u>) = u>rij(u))/c 
is the associated wave number. 

The Hamiltonian (||) yields the state at the output of 
the crystal (for vacuum-input signal and idler) via, 
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df 5/(010), 
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for small values of the coupling constant x*" 2 "*- As done in 
[pi, we shall assume that x' 2 ' is independent of frequency 
over the pump bandwidth, even though this assumption 
may not be satisfied in some ultrafast applications. For a 
strong coherent pump pulse and in the absence of pump 
depletion, we may replace the pump field operator in 
Eq. (B) with, 
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where £ p (u>) is a classical complex amplitude. Because 
we are interested in the fields far from the crystal, we may 
expand the integration limits in Eq. ([5]) to run from — oo 
to +oo. Thus, the t' integration produces an impulse, 
6(ujp — lj s — LUi), that expresses energy conservation at 
the photon level. 

The biphoton state that we arc seeking is the non- 
vacuum part of Eq. (^|). Under the preceding assump- 
tions, this is given by, 
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where = a^(w)|0) is a single-photon state, 



a(uj s ,LJi) = — UJsUJl — -£ p (Lj s +u!i), (8) 
n s (u> s )rii(u)i) 



is determined by the pump spectrum, and 
, sin(Afc(o; s ,^)£/2) 

is the phase-matching function, with Ak{uj s ,uji) = 
k p (uj s + uji) ~ k s (uj s ) ~ ki(t0i). 

To obtain maximal entanglement from the biphoton 
state (0) we need to collapse the double integral over fre- 
quency into a single integral. For the customary twin- 
beam state |TB), this is accomplished by using a cw 
pump of frequency lo p to force a (u> s , u>i) oc 5(u> s +u>i — u) p ) 
in Eq. (Q). We then obtain a TB state (1) with spectral 
function <f>{oS) oc &l(w p /2 — lo,lo p /2 + ui). Because this 
makes the common signal/idler fluorescence bandwidth, 
ilf, inversely proportional to L, it follows that short crys- 
tals are better suited to generating broadband TB states. 
For DB states, however, we will see that long crystals do 
not prevent broadband entanglement generation. 

Continuous- wave operation is not the only way to ob- 
tain a maximally-entangled state from (0). We can 
also eliminate one of the frequency integrals by forc- 
ing to approach a delta function. The prop- 
erty limi^ OQ [sin(a;L)/x] = ttS(x) allows us to write 
<J>l(oj s , u>i) — 2ttS (Ak(u) s ,Wi)) for an infinitely long crys- 
tal. (In practice the nonlinear crystal will always have a 
finite length L, but we will see that a high degree of en- 
tanglement can be obtained using practical values of L.) 
To force <J>l(w s , u>i) oc 5(ui s — uji) in the long-crystal limit, 
we need to ensure that Afc(cj s , cjj) = if and only if lo s = 
uii, for lo s +u>i ranging over the full pump bandwidth, £L,. 
Equation (Q) then reduces to the DB state of Eq. (g), 
with spectral function <fi(u>) = a{u> p /2 + u,lo p /2 + oS), 
where lj p is the pump beam's center frequency. Note 
that <fi(u>) depends only on the pump spectrum and the 
refractive indexes of the nonlinear crystal, as can be seen 
from Eq. (||), and that its bandwidth is Q p /2. Moreover, 
the symmetry of the phase-matching function forces 
the signal and idler fluorescence spectra to be identical, 
something that is not generally true in ultrafast type-II 
downconversion ||. 

Is it possible to satisfy the condition Ak(u> s , u>i) = 
only for to s = Ui over the full pump bandwidth? What 
does this condition correspond to physically? By using 
the first-order Taylor expansions of k s and ki around 
cjp/2 and of k p around co p , we find that, 



n p (u p ) 
k' p {u p ) 



n s (uj p /2) + nj(uj p /2) 
2 : 
^K/2) + ^K/2) 



(10) 

(11) 



ensure that Ak(oj p /2 + lj,uj p /2 + uj) = for \lj\ < O p /2. 
In physical terms, the extended phase-matching condi- 
tion given by ( |l(]| ) and (11) assert that the index of re- 
fraction and the inverse group velocity seen by the pump 
must equal the averages of those seen by the signal and 
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idler. Equation ( |To| ) is the customary phase-matching 
condition required for the generation of |TB) at fre- 
quency degeneracy. Equation ([ll]) is equivalent to the 
"group velocity matching" condition introduced in ||. 
It turns out, however, that Eqs. ([!(]) and ( JPI| ) are not 
sufficient for DB state generation. Because Ak(uj s ,uji) 
must vanish only for oj s = u>i, we must also require that 
k' s (ujp/2) ^ k' i (uj p /2). This requirement excludes type-I 
crystals, for which k s (oj) = ki(uj). Thus, in all that fol- 
lows we will presume type-II operation. We will discuss 
later the validity of truncating the Taylor series at the 
n = 1 terms. 

The states |DB) and |TB) are duals in the follow- 
ing sense. The former is a biphoton whose constituent 
photons are coincident in frequency, and the latter is a 
biphoton whose constituent photons are time-coincident. 
We now show that coincidence counting using Hong- 
Ou-Mandel (HOM) and Mach-Zehnder (MZ) interfer- 
ometers, as sketched in Fig. [j], can provide experimen- 
tal quantum-interference signatures that distinguish be- 
tween the TB and DB biphoton states. 



± cos[(wi ± w 2 )t] J . 



(13) 
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FIG. 1: Quantum interference experiments to distinguish 
jDB) from |TB). The left panel shows the HOM interfer- 
ometer and the right shows the MZ interferometer. In both 
cases, coincidence counts are measured as the relative delay, 
r, between the interferometer's arms is varied by moving the 
beam splitter that is nearest to the detectors. The half- wave 
plate (HWP) rotates the idler polarization to match that of 
the signal because type-II downconversion is assumed. 

To understand the outcome of the two experiments 
shown in Fig. |l|, it is useful to start from the general state 
| \&) of the form (Q). The coincidence rate for detection 
intervals that are long compared to the reciprocal fluo- 
rescence bandwidth is given by (see for example [[j], 0), 
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where a± and a 2 are the photon annihilation operators 
at the two detectors and |\&) is the biphoton state of the 
source. Assume that the product a(u>i, W2)^l(wi, ^2) is 
symmetric in uji,^, as is the case for both |TB) and 
|DB). By applying the beam-splitter transformations on 
the operators a s and ai, we find that, 

P±{t) oc J^- /~2^ H w i> w 2) $l(^i,^2)| 2 



In (|13|) the minus signs apply to the HOM interferometer, 
and the plus signs apply to the MZ interferometer. 

For the TB state, we set \a(u>i , 0J2) | 2 oc S(ui p — uji — 1^2) 
and |$ L (wi,w 2 )| 2 oc |<£[(w 2 - ^i)/2]| 2 . Ap- 
proximating the fluorescence spectrum by 
\(/)(uj)\ 2 = sm 2 (2iruj/n f )/[2iruj/(n f L)} 2 , where 
fl f = 4ir/(L\k' s (Lu p /2) - fej(w p /2)|), Eq. © then 
yields the familiar triangular-shaped HOM coincidence 
dip of width Air/Q. f centered at r = |[ g|. The 
|TB) coincidences require that one photon exits from 
each output port of the beam splitter. At zero relative 
delay, the two quantum trajectories that give rise to 
such coincidences destructively interfere, leading to a 
coincidence null The width of this dip is ~l/f2/, 
because signal and idler wave packets separated by many 
reciprocal fluorescence bandwidths are distinguishable 
and hence do not interfere. When the TB-state coinci- 
dence rate is evaluated for the MZ interferometer, we 
obtain P+{t) cx 1 + cos(o;pr), i.e., sinusoidal fringes at 
the pump frequency. These fringes have infinite extent 
because a perfect cw pump has infinite coherence time. 

Now suppose that the input state in Fig. [I] is |DB), 
i.e., let |$l(wi,w 2 )| 2 oc 5(uji - u 2 ) and |a(wi,w 2 )| 2 = 
|0[(wi + UJ2 — w p)/2]| 2 in Eq. (p"3|). In this case the fre- 
quency coincidence between the signal and idler photons 
eliminates any delay dependence in the HOM configu- 
ration, reducing Eq. (|T^) to P-(t) = 0. In fact, the 
wavefunctions for the two photons extend to all times 
and cannot be separated: the quantum trajectories that 
give rise to coincidences destructively interfere for any 
delay r. For the Mach-Zehnder arrangement, the DB 
state gives P+(t) oc 1 + exp(— 17 2 r 2 /4) cos(w p r), under 
the assumption of Gaussian pump spectrum |0(w)| 2 oc 
cxp[— 4w 2 /f2 2 ]. Notice that in this case P+ again ex- 
hibits pump-frequency interference fringes, but now the 
interference pattern has width A/Vt p , i.e., roughly equal 
to the duration of a transform-limited pump pulse. Sim- 
ilar interference patterns have been previously analyzed 
m|. 

Both the HOM and the MZ interferometers distinguish 
between the states |DB) and |TB). However, because DB 
state generation requires infinite crystal length — whereas 
TB state generation uses a finite-length crystal — it be- 
hooves us to study what happens in the finite-L regime 
when we use a pulsed pump in conjunction with our ex- 
tended phase-matching conditions. The biphoton state, 
|DBl), that this system generates is entangled in fre- 
quency, but not maximally so, i.e., measuring the fre- 
quency of the signal photon does not exactly determine 
the frequency of the idler photon. When |DBi) is mea- 
sured with an HOM interferometer, the resulting coin- 
cidence null is no longer of unlimited extent. Indeed, 
the width of the coincidence dip for |DBi) is identical 
to that for |TB). Thus the HOM interferometer cannot 
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distinguish between these two states. The MZ interfer- 
ometer, however, does distinguish between |DBl) and 
|TB), as the width of the former's fringe pattern is set 
by the pump bandwidth and hence independent of crystal 
length. 

More insight into the complementary behavior of |TB) 
and |DB) can be gained by examining their time do- 
main structures. Both of these biphoton states arise from 
the coherent superposition of spatially-localized, instan- 
taneous signal/idler pair creations occurring throughout 
the length of the nonlinear crystal. HOM and MZ in- 
terferometers use integrating photodetectors, but reveal 
temporal aspects of |TB) and |DB) via quantum inter- 
ference. Suppose, however, that we use an ultrafast pho- 
todetector to measure the arrival time of the signal pho- 
ton. This measurement specifies a definite location, along 
the crystal, at which the detected signal photon was cre- 
ated, and implies rather different temporal statistics for 
its associated idler photon depending on whether the 
biphoton was |TB) or |DB). Because the TB state is 
produced by a cw pump, its component photons may be 
created at any time. However, once its signal photon has 
been detected at time T s , the accompanying idler photon 
must be at Tj, where \T S — Ti| < 47r/f2/ for our type-II 
system. The individual photons in the DB state also may 
be created at any time, even though this biphoton is gen- 
erated by a pulsed pump. Here, the timing uncertainty 
is really uncertainty in the location, within the infinitely- 
long crystal, at which the photon pair is generated. Once 
again, detection of a signal photon at time T s provides 
location information which strongly constrains the ar- 
rival time for the idler photon. In particular, the ex- 
tended phase-matching conditions that produce the DB 
state under pulsed pumping force (T s + Ti)/2 to have a 
mean value at a fixed offset — set by dispersion — from the 
peak of the classical Gaussian pump pulse. 

It turns out to be difficult to find a crystal satisfy- 
ing the two conditions ([Io|) and (|IT|). Thus, we will 
enforce (|1Q) via quasi-phase-matching in a periodically- 
poled x material Efl], i.e., one for which the addi- 
tion of an artificial grating results in a spatially-varying 
nonlinear coefficient, X ( z ) — exp(i27rz/A), along 
the propagation axis. By choosing the grating period 
A to cancel the zeroth-order term in the Ak(u s ,uJi) ex- 
pansion, we can replace Eq. ([l(]) with the new con- 
dition n p (ujp) — [n s (ujp/2) + rii(ujp/2)]/2 — 2irc/(Acu p ). 
This, together with Eq. ([ft]) is satisfied by PPKTP at 
a pump wavelength of 790 nm with a grating period of 
47.7/zm when propagation is along the crystal's X axis, 



the pump and idler are ^-polarized, and the signal is Z- 
polarized. It still remains for us to examine the validity 
conditions for the L —* oo approximation to the phase- 
matching function $j,(u; s , caj). These can be shown to be 
2ir/-yn p < L < 8tt//x^, where 7 = \k p (uj p ) - k' s (cj p /2)\, 
and /x is the maximum-magnitude eigenvalue of the Hes- 
sian matrix associated with the 2-D Taylor series expan- 
sion of Ak(u) s ,tUi). Physically, the lower limit on crystal 
length is set by our need to be in the long-i regime, and 
the upper limit is set by the second-order terms in the 
Taylor expansion. For our PPKTP example, we have 
that 7 w 1.4 x 10~ 4 ps/^m and fi w 3.6 x 10 _7 ps 2 //um. 
With a 170 fs (fl p /2ir = 3THz) transform-limited pump 
pulse the preceding crystal-length restrictions reduce to 
0.23 cm <C L <C 19.7 cm, so that a 2-cm-long crystal will 
suffice. Finally, we note that polarization-entangled DB 
states can be created by paralleling the procedure in |Tl|l 
for the creation of polarization-entangled TB states. 
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